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Abstract 

We present a survey on several mass transportation problems, in which a given 
mass dynamically moves from an initial configuration to a final one. The approach we 
consider is the one introduced by Benamou and Brenier in [5j, where a suitable cost 
functional F(p, v), depending on the density p and on the velocity v (which fulfill the 
continuity equation), has to be minimized. Acting on the functional F various forms 
of mass transportation problems can be modeled, as for instance those presenting con- 
gestion effects, occurring in traffic simulations and in crowd motions, or concentration 
effects, which give rise to branched structures. 
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1 Introduction 

Mass transportation theory goes back to Gaspard Monge: in 1781 he proposed in [24J a model 
to describe the work necessary to move a mass distribution /ii = p\ dx into a final destination 
A*2 = P2 dx, given the unitary transportation cost function c(x, y) which measures the work 
to move a unit mass from x to y. The goal is to find a so-called optimal transportation map 
T which moves Hi into fi2, i-e. such that 

Hi{E) = /ii(T _1 (_E)) for every measurable set E, 

with minimal total transportation cost 

J c(x, T(x)) dpi. 

The measures /ii and \i2 have equal mass (normalized to one for simplicity) and are called 
marginals; the operator 

T*iM(B) = fi(T-\B)) 
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is called push forward operator. The optimization problem then becomes 



min | J c(x, T(x)) dp,\ : T* p,\ = /i 2 |. 



The natural framework for this kind of problems is the one where X is a metric space and 
/ii,/i2 are probabilities on X; however, the existence of an optimal transport map is a very 
delicate question, even in the classical Monge case, where X is the Euclidean space M. d and 
c(x,y) = \x — y\ (see for instance [2S1 SI [201 US, [26J ) . Thus in 1942 Kantorovich proposed 
in [22] a relaxed formulation of the Monge transport problem: the goal is now to find a 
probability on the product space, which minimizes the relaxed transportation cost 

C(/i , yUi) = J c(x, y) 7(dx, dy) (1.1) 

over all admissible probabilities 7 on X x X, where admissibility means that the projections 
7r*7 and 71"* 7 coincide with the marginals fii and ^2 respectively. The Kantorovich problem 
then reads 



mm 



c(x, y) 7(cfe, dy) : vr # 7 = fij for j = 1, 2 



The cases c(x, y) = \x — y\ p with p > 1 have been particularly studied, and the cost C(/i , /ii 



in (1.1) provides, through the relation 

Wp(/j,o, fii) = (C(/io,/ii)) 1/p , 

the so-called Wasserstein distance w p which metrizes the weak* convergence on the space 
of probabilities V(Q). A very wide literature on the subject is available; we mention for 
instance the books [21 [2Z1 12H] where one can find a complete list of references. 

A dynamical model of mass transportation consists in finding, given two probabilities po 
(the initial configuration) and p\ (the final destination), a curve pt of probabilities, with t G 
[0, 1], joining po to p\ and minimizing some suitable functional. This functional should take 
into account the properties of the particular transport under consideration, like congestion 
or concentration phenomena. 

Several models have been proposed to describe the various effects that may occur during 
a transportation path; in the next sections we will describe some of them. We want to 
stress the fact that in the congested situations (like for instance the crowd motion in case of 
panic) people move looking for paths with a low mass density, because they allow a higher 
velocity; we will see that this feature can be well described by a convex cost functional. On 
the contrary, in the cases where a mass concentration is favoured (as for instance in public 
transportation networks), people look for paths with a high mass density, which is related 
to some concavity in the cost. 



2 The path functionals approach 

In this section we consider a very general way to define minimizing paths for functionals 
defined on curves in a metric space (X, d). The interesting case for dynamical models in 
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mass transportation is when X coincides with the class of all probabilities V(Q) endowed 
with the Wasserstein distance w p . 

Let (X, d) be a metric space; for simplicity we assume that all closed bounded subsets of 
X are compact. For every Lipschitz curve 7 : [0, 1] — > X we define the metric derivative of 
7 at the point t as 

By Rademacher Theorem it can be seen (see [3]) that for a Lipschitz curve 7 the metric 
derivative exists for a.e. t G [0, 1] and that 



d(7(*),700) < / |Y|x(r)dT Va,f e [0,1]. 
We are concerned with variational problems for functionals of the type 

J(i)= I J{i{t))\i\ x {t)dt (2.1) 



where 7 : [0, 1] — > X varies among all Lipschitz curves with fixed endpoints 7(0) = x and 
7(1) = Xx. 

Theorem 2.1. Let X be a metric space such that all closed bounded subsets of X are 
compact, let J : X — )■ [0, +00] be a lower semicontinuous functional on X, and let x , X\ G X 
be fixed. Assume also that 

(HI) the functional J is finite on at least a Lipschitz curve 70 joining xq to x±; 
(H2) the functional J is bounded from below by a constant a > 0. 
Then the minimum problem 

min|j7"(7) : 7 Lipschitz, 7(0) = x , 7(1)= 07 j (2.2) 

admits a solution. 



The assumption of Theorem 2.1 on the metric space (X, d) can be slightly relaxed and 
the following result holds. 

Theorem 2.2. Let (X, d) be a metric space and let d' another distance on X such that: 
(Kl ) d! < d; 

(K2) all d-bounded sets in X are relatively compact with respect to d' ; 

(K3) the mapping d : X x X — )■ 1R + is lower semicontinuous with respect to the distance 
d! x d'. 



Let J be the functional defined in (2.1), where J : X — > [0, +00] is assumed lower semi- 
continuous with respect to d' . Then, under the assumptions (HI) and (H2) above on J, for 
every Xq,X\ G X the minimum problem (2.2) admits a solution. 
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Remark 2.3. The coercivity assumption (H2) of Theorem 2.1 above can be weakened requir- 
ing only that 

( inf J) dr = +00 

V B r (x) / 

for some (hence for all) x G X. Assumptions (HI) and (H2) above can still be weakened by 
simply requiring that there exists an admissible Lipschitz curve 70 such that 

"+00 



J(lo) < [ ( inf j) dr. 

J \B r (x ) J 



We refer to [9] for the proofs of the results above; we apply here these results to the 
case when X is a Wasserstein space of probabilities. More precisely, we consider a compact 
metric space Q equipped with a distance function c and a positive finite non-atomic Borel 
measure m. We consider the q- Wasserstein metric space W q (Q) of all probability measures 
juonfi, equipped with the p- Wasserstein distance (with q > 1) 



WqQi,!, H2) = inf ( / c(x, y) q X(dx, dy) 



1/9 



where the infimum is taken on all transport plans A between /i x and fi2, that is on all 
probability measures A on Q x Q whose marginals nf A and 7rf\ coincide with \i\ and /i 2 
respectively. 

We consider functions J : W 9 (f2) — >■ [0, +00] to be used in (2.1). Functionals of this 
form have been studied by Bouchitte and Buttazzo in a series of papers (see [TJ |8]) and 
it is shown that, under the weak* lower semi continuity and a locality property, they can be 
represented in the integral form 



J( M ) 

where 



m is a nonnegative nonatomic finite measure on Q; 

f : Q x R — > [0, +00} is an integrand with f(x, •) convex and lower semi continuous; 
dfi/dm is Radon-Nikodym derivative of \x with respect to m; 

/i s is the singular part of \i with respect to m according to the Radon-Nikodym de- 
composition theorem; 

/°° is the recession function of /, defined by (the limit is independent of the choice of 
s in the domain of f(x, ■))): 

r (*,.)= iim /(x - s ! + is) ; 

is the set of atoms of /i, i.e. the points iffi such that fi(x) := > 0; 
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• # is the counting measure; 

• g : ft x R — > [0, +00] is an integrand with g(x, •) subadditive and lower semicontinuous, 
such that g(x,0) = 0, and fulfilling the compatibility condition 

, . g(x,st) rno , . 

g (x,s)= sup — = / {x,s). 

t>0 t 

Example 2.4. Taking f(s) = \s\ p with p > 1 and = +00 (with g(0) = 0) we have the 
Lebesgue type functionals 



JQi) 



w| p dm if /i = u dm with it G L p (f2, m) 
00 otherwise 



whose domain is L p (Q,m). On the other hand, taking f(s) = +00 and g(s) = \s\ a with 
a < 1, we have the Dirac type functionals 



Jin) 



EKI"= if = E t o»ix» is a discrete measure 



k 

-00 otherwise 



whose domain consists of discrete measures. Finally, taking f(s) = \s\ p with p > 1 and 
g(s) = with a < 1 we have the Mumford-Shah type functionals 

J \u\ p dm + J \n(x)\ a d#(x) if fj, — u dm + J2k a ^x k 
I +00 otherwise 

whose domain consists of measures with no Cantor part. 

Putting together the characterization of weak* lower semicontinuous functionals J(n) 
above, and the abstract existence Theorem 2.1 we have the following result (see [9] for the 
proof). 

Theorem 2.5. Suppose that f(s) > for s > and that g(l) > 0. Then we have J > c > 0, 
so that the coercivity condition (H2) of Theorem 2.1 is fulfilled. Therefore, the minimum 
problem (2.2) admits a solution, provided that there exists at least a Lipschitz curve 7 : 
[0, 1] — > W q (Q), with given starting and ending points, with finite cost. 

It is interesting to study the situations when, given two probabilities fiQ and fii on an 
Euclidean domain Q C lR d , an optimal path joining them actually exits. We analyze in a 
more detailed way the two cases below. 

• f(s) = \s\ p with p > 1 and g(s) = +00 (with g(0) = 0), which can be used to describe 
the congestion phenomena; the functional (2.1) then becomes 

-1 



J P (l) = J o ( JJl(t)\ p dx)W\ Wq (t)dt 



where W q is the g-Wasserstein space and the expression J n \ j{t)\ p dx has to be intended 
as +00 when the measure ^(t) is singular with respect to the Lebesgue measure dx. 
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f(s) = +00 and g(s) = \s\ a with a < 1, which on the contrary describes the concen- 



tration phenomena. In this case the functional (2.1) takes the form 



Ja{l) = J 1 (Jjj(t)\ a d#)\j'\ Wv (t)dt 



In the first case the domain of the functional J is L P (Q) and the question is answered by 
the following result (see |9J). 

Theorem 2.6. If Ho = uo£ d o,nd n\ = u\L d , with «o, «i £ L p iS^), then /io and fi\ can always 
be joined by a finite energy path, hence by an optimal energy path. 

In addition, if p < 1 + then any pair of probabilities fio, A*i can be joined by a finite 
energy path, hence by an optimal energy path. 

Finally, if p > 1 + 1/d, every nonconstant Lipschitz path ^(t) starting from a Dirac mass is 
such that Jpi^f) = +00. Therefore, if p > 1 + 1/d no measure can be joined to a Dirac mass 
by a finite energy path. 

The second case, occurring when concentration phenomena are present, is somehow sym- 
metric; in this case the domain of the functional J is made only of discrete measures and 
the existence of optimal paths is solved by the following result (see [§]). 

Theorem 2.7. // fiQ and fii are finite sums of Dirac masses, i.e. /io = cii5 Xi and 
fii = ^2,jbj5 yj , then fiQ and /ii can always be joined by a finite energy path, hence by an 
optimal energy path. 

In addition, if a > 1 — 1/d, then any pair of probabilities /io,/ii can be joined by a finite 
energy path, hence by an optimal energy path. 

Finally, if a < 1 — 1/d, every nonconstant Lipschitz path 7(t) starting from the Lebesgue 
measure is such that Jai^f) = +00. Therefore, if a < 1 — 1/d no measure can be joined to 
the Lebesgue measure by a finite energy path. 

The example below shows a simple case in which an explicit computation of the optimal 
concentration path can be made. 

Example 2.8. In the Euclidean plane M. 2 consider the origin O, the points A = (l,a) and 
B = (1,— a) with a > 0, and take /xq = So, A*i = \Sa + \Sb- Using the concentration 



transport functional of Theorem 2.7 with < a < 1, an easy calculation shows that the 
optimal path j(t) is given by 



7(*) 



S(t,o) for < t < r a 



\S x {t) + \S y {t) for r a < t < 1 
where r a = (l - a(4 x - a - l)~ 1/2 ) + and 

x(t) = (t,a - y(t) = _ T " )- 



Figure [T] shows the plot of the entire path in the case a = 1 for some values of the parameter 
a. Notice that in this = for a > 1/2, while for a = we recover the Steiner 

problem. 
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Figure 1: The optimal concentration paths for a = 0.5, a = 0.25, a = 0. 

3 The Eulerian approach 

A different dynamical approach to mass transportation problems was proposed by Benamou 
and Brenier in [5] (see also [15]). It consists in looking to pairs (pt,v t ) where p t represents 
the mass distribution at time t and v t the related velocity field; from the mass conservation 
the pair (p, v) has to satisfy the so-called continuity equation 

d t p + div x (pv) = 

Among all pairs (p, v) that satisfy the continuity equation above, the one providing the 
mass transportation, from an initial configuration po to a final one pi, is obtained by the 
minimization of a suitable functional J-(p, v) that in [5] is taken equal to the kinetic energy: 



F(p,v) = J (J \v\ 2 dp^j dt. 



In this way the probability po is dynamically transported on the probability p\ following the 
geodesic path on the 2-Wasserstein metric space W2. 

To be more precise, by Theorem 8.3.1 of (2] we have that for every absolutely continuous 
curve p t in the p-Wasserstein metric space W P (Q) (with p > 1) there exists a map q from 
[0, 1] into the space of vector valued measures, such that q t <C p t (hence q t = v t pt, being v 
the velocity vector) which represents the flux q = pv and satisfies 

d t p + div x q = and ||^t||iP( pt ) = \p' t \w p for a.e. t G [0, 1]. (3.1) 

The kinetic energy is replaced in this case by the action functional 

FM= f( f P- P dpt)dt. (3.2) 



dp t 

In the degenerate case p = 1 a little more care is needed, since the absolute continuity q t pt 
is no more guaranteed, and the L 1 -norm has to be replaced by the mass of the measure q t 
(see HI). 
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Note that, using the variables p and q instead of p and v , provides the convexity of the 



functional T(p, q) in (3.2). 



The precise meaning of the continuity equation 

d t p + diyj. q = in [0, 1] x Q 
q ■ v = on [0, 1] x dQ, 

has to be given in the sense of distributions, that is 



d t (j)(t,x) dp t (x) 



D x (j)(x,t) ■ dq t {x) 



dt = 



for every smooth function with 0(0, x) = 0(1, x) = 0. 

The general dynamical formulation of a mass transportation problems, following this 
Eulerian formulation, then becomes 



min{j r (p,g) : d t p + div x q = 0, p(0, •) = p , p(l, •) = Pi}. 



(3-3) 



In the minimization problem above the continuity equation provides a linear constraint, and 
the existence of an optimal dynamical path p t easily follows by the direct methods of the 
calculus of variations. In the theorem below (see [16]) we denote by Q the time-space domain 
[0, l]xfi C IR 1+d , with outer normal versor n, and by a the measures of the form (p, q), which 
belong to the space Aib{Q, M 1+d ) of IR 1+d -valued measures defined on Q. The minimization 



problem (3.3) can be then written in the form 
min {J-{cr) : 



div a = f in Q, a ■ n = on dQ } , 



(3.4) 



where the scalar measure / = 5\(t) <8> pi(x) — 5o(t) ® po(x) takes into account the initial-final 
configurations. 

Theorem 3.1. Let J 7 : Aib{Q, M. 1+d ) — > [0, +oo] be lower semicontinuous for the weak* 
convergence, and assume that the coercivity condition 



Ho) > c \°\ 



1 

c 



W e M b (Q, 



(3.5) 



holds for a suitable constant C > 0, where \o~\ denotes the total variation of a on Q. As- 
sume also that J-{(Tq) < +oo for at least one measure a satisfying the continuity equation 
— diver = f inQ, with the boundary conditions o~-n = on dQ. Then the minimum problem 



(3.4) admits a solution. Moreover if J- is strictly convex, this solution is unique. 



In the case when T is convex, problem (3.4) also admits the dual formulation 

min { Tip) : — diva = / in Q, a-n = on 9Q} 
sup{ / V df-T(D V ) : v eC\Q)}, 



(T*oA)*(f) 



where A : C{Q) — > C(Q,M> 1+d ) denotes the operator with domain C X {Q) given by 

A{ip)=D^> V^GC^Q). 
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The dual formula above holds if J 7 * is continuous at least at a point of the image of A. The 
primal-dual optimality condition then reads as 

J Difopt ■ dcr opt = ^{(Topt) + jr*(D(popt), (3.6) 

provided an optimal solution ip opt of the dual problem exists. The point is that, in general, 
the maximizers ip op t of the dual problem are not in C 1 (Q) and a relaxation procedure is 
necessary to make the primal-dual optimality condition meaningful. We do not deal with 
this rather delicate question, and we refer the interested reader to [To] . 

Several variants of mass transportation problems have been studied by considering in 



(3.3) various convex functions of the pair (p,q). 

• Dolbeault, Nazaret and Savare considered in UHl functionals of the form 



Hp, q) = J ( J $G°> ?) dm ) dt > 



where m is a given reference measure on M. d , p and q are identified through their 
densities with respect to m, and 



Note that, if the function h is concave (for example h(p) = p^, with e [0, 1]), the 
functional J 7 turns out to be convex as well. Functionals of this kind are motivated to 
provide efficient models to study diffusion PDEs of the nonlinear mobility type 

dtp + div x (h(p)v) = 0, 

which can be interpreted as gradient flows of a given functional with respect to a new 
family of distances generalizing the ones of Wasserstein type. 



Buttazzo, Jimenez and Oudet considered in [16] a model to describe the behaviour of 
a crowd under some panic effects. The dynamical model is as above, with 

jF{p, q) = J ^ J — + cp 2 dx^j dt, 

which is a convex functional defined on Lebesgue integrable densities (p,q). 

Maury, Roudneff-Chupin and Santambrogio considered in [23] the problem of an effi- 
cient emergency evacuation of a crowd; the target configuration p\ is not prescribed, 
being replaced by an integral functional cost which takes into account the goal of the 
crowd. Given a room Q with an exit T out the model consists in a gradient flow evolu- 
tion, in the 2- Wasserstein metric space W2(fi), from an initial given density po- The 
functional governing the model is 



Hp) 



d(x) dp if p G K 

n 

+oo otherwise 
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where d(x) denotes the distance of the point x G Vt from T out and 



K = {pe V(Q), p = pa(x)dx + p out , p n < 1, supp(p OMi ) C T out }. 

The situation is more delicate if we want to provide model mass trasportation models 
with concentration effects through an Eulerian formulation involving the continuity equation. 
This because we have to consider in this case concave functionals on the space of measures, 
taking into account the concentration effects. This happens when the moving mass has the 
interest to travel together as much as possible, in order to save part of the cost; for instance 
this occurs in the transportation of signals along telephone cables, as discovered by Gilbert 
who in [21] formulated a mathematical model for it. 

The starting point is, as above, to consider the admissible class T> of all pairs (p, q) with 
p G C([0, 1);V(Q)) and q G ^([0, 1]; M(Q; R d )) satisfying the continuity equation 

d t p + divr q = in [0, 1] x Q 
q ■ v = on [0, 1] x <9f2, 

and the class D(po,pi) of all pairs in V with initial and final data p(0, •) = po, p(l, •) = p\- 
Our goal is to minimize on "D(p ,pi) an integral functional of the form 

HP' 0)= F (Pt> It) dt. 
Jo 

In [T2] we presented some natural choices for the function F(p t , qt) and we showed that the 
integral functional J 7 above is both lower semicontinuous and coercive with respect to a 
suitable convergence on (p, q), which directly provides the existence of an optimal dynamical 
path. 

Denoting by G a (0 < a < 1) the functional defined on measures 



G a (X) 



G a {\v\ l / a ■ p) iiq = vp, 
-oo if q is not absolutely continuous w.r.t. p. 



I \\({x})\°<d#(x) = Y,\M a ii\ = J2 ien kS Xi 

+oo if A is not purely atomic, 

we define 

In this way our functional T becomes 

HP,Q) = I \ I \vt{x)\ Pt {{x}) a d#{x)] dt= [ \J2MpI] dt, (p,q)eV, 
Jo l Jn J Jo 1 teN J 

and the dynamical model for branched transport we consider is 

B a (po,Px) ■■= min {T(p,q) : (p, g) G £>(p , Pi)}. (3.7) 
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Remark 3.2. We point out that the weak* convergence of the pairs (p,q) is too weak for 



our purposes; indeed it does not directly imply the lower semicontinuity in (3.7), since the 
functional F is not jointly convex in the pair (p,q). On the other hand, if {p n ,q n ) is a 
sequence in V, and we assume that 

(j>t,q?)^(j>t,qt), for^-a-e. te[0,l], 

then a simple application of Fatou's Lemma gives the desired semicontinuity property of 
F. In fact, as a consequence of the semicontinuity of G a and of the convexity of the map 
(x,y) i — y \x\ p /y p ~ l , we obtain that F is a lower semicontinuous functional on measures. 



In order to prove the existence of minimizers for problem (3.7) through the direct methods 
of the calculus of variations involving semicontinuity and coercivity results, we introduce a 
convergence which is stronger than the weak* convergence of measures on [0, 1] x Q, but 
weaker than weak* convergence for every fixed time t. 

Definition 3.3. We say that a sequence (p n , q n ) in T> r-converges to (p, q) if(p n , q n ) — ^ (p, q) 
in the weak* sense of measures and in addition 

sup (F(p"g t n ) : n6N, t G [0, 1]} < +oo. 

Note that, due to the fact that the functional F is 1-homogeneous in the velocity, its value 
does not change by a reparametrizations in time. By reparametrization we mean replacing 
a pair (p,q) with a new pair (p,q) of the form p t = p v (t), qt = ^p'{t)qtp{t)- This equivalently 
means that q is the image measure of q through the inverse of the map (t,x) i-» (ip(t),x). 
The following results are obtained in [12J. 



Theorem 3.4. Let (p n , q n ) be a sequence in V such that F(p n , q n ) < C for a suitable 
constant C. Then up to a time reparametrization, (p n ,q n ) is r-compact. 

Theorem 3.5. Let (p n ,q n ) be a sequence in T> r-converging to (p,q). Then 

F(p,q) < liminf F{p n ,q n ). 

n— >oo 

As a consequence we obtain the following existence result. 



Theorem 3.6. For every po, Pi G F{Vt), the minimization problem (3.7) admits a solution 

( P ,q)ev. 

Remark 3.7. It has to be noticed that, similarly to what happens in the path functional ap- 
proach of Section [2j for some choices of the data po, p\ and of the exponent a, the functional 
F could be constantly +oo on every admissible path (p,q) G V(p , pi) joining p to p\. In 
[T2] the equivalence of the Eulerian model above with other variational models for branched 
transportation has been proven. For these models finiteness of the minima has been widely 
investigated, so we can infer for instance that if a > 1 — 1/d then every pair po and pi can 
be joined by a path of finite energy. On the other hand, if a < 1 — 1/d, p = S xo and pi 
is absolutely continuous with respect to the Lebesgue measure C d , then there are no finite 
energy paths connecting them. 

Some configurations of branched transportation paths have been computed numerically 
by E. Oudet, providing the outputs of Figures [2] and [3] (for more examples we refer to the 



web site http : / /www . lama . univ-savo ie . f r/~oudet| ). 
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Figure 2: Optimal transport path of a point in 4 points: a = 0.6, 0.75, 0.85, 0.95. 



4 Other models 

In addition to the evolution models presented in the previous sections other approaches 
are possible. We present here a model which is under investigation in |17j . It consists in 
assuming that the mass is composed by many particles, each one moving under the action 
of a potential which takes into account the mutual interaction among them. 

If Xi(t) is the position of the z'-th particle at time t the Lagrangian function governing 
the system is 

L(X,X>) = ^Y,\ X 'i\ 2 + W^ V{Xi > Xj) 

i ijij 

where N is the number of particles, i,j — l,...,N, and V denotes the potential describing 
the mutual interaction among particles. 

When the number N of particles goes to infinity, instead of describing the motion of 
every single particle, it is interesting to consider the Eulerian formulation which consists in 
describing the movement of the mass density p(t,x), which satisfies the continuity equation 

d t p + divq = 0. (4.2) 

The pair (p, q) will then be seen as the minimizer of a suitable functional, obtained by passing 



to the limit as iV — > +oo the Lagrangian functional (4.1). 

The potential V will of course enter in the expression of the limit functional. The term 
2^ l^il 2 produces the kinetic energy J y- dx, while the potential part V(Xh Xj) 

— P 12 ' 




Figure 3: Optimal transport path of a point in a circle: a = 0.6, 0.75, 0.85, 0.95. 

gives rise to the term ff V(x,y)p(x)p(y) dxdy. Summarizing, we end up with the energy 
functional ^ 

e(p, q) = J {f^2p~ dx + ff v ^ x ' y)p( x )p(y) dx d y) dt 



that has to be minimized among all pairs (p,q) which satisfy the continuity equation (4.2) 
together with initial and final conditions p(0, •) = po an d p(l, •) = pi- 

The energy E(p, q) above has to be better detailed when p(t, •) is a singular measure, 
situation that may happen with some potential V. In the singular case the expression of 
E(p, q) becomes 

E( P ,q) = £ {If \jp\ dp+ S I v ^y)p^ dx )p^y)) dt ( 4 - 3 ) 

Under some mild conditions on the potential V it is possible to obtain an existence result 
for an optimal path p(t,x). Its behaviour, in terms of congestion or concentration effects, 
heavily depends on the form of the potential V. Clearly, a potential of repulsive type will 
produce a congestion effect, while a potential of attractive type will produce concentration. 

Example 4.1. In the two-dimensional case assume that the initial and final configurations 
are 

1 r 1 r 

Po = o A , Pi = + 7^, 
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where A = (0,0), B = (1,1), C = (1,-1). In the case of an attractive potential the 
evolution consists of the motion of two particles that travel together up to a certain point 
and then move symmetrically to reach their final destinations, respectively the points B and 
C. If (x(t),y(t)) is the path of the upper particle, we have x(t) = t and y{t) minimizes the 
functional 



f 1 (\yf + V(2y))dt. 
Jo 



Then y(t) — for t <t , where 



to = 1 



{V(2y)y 1/2 dy, 



while 



y " = W(2y) for t > t . 
In the case V(y) = C\y\ a with a > we have 

2 (2-a)/2 

to = if a > 2, t = 1 - — t= if a > 2. 



If < t < 1 we find 



(2-a) v / C 



//(/) '(VYV2- n)(/ / n )) 2:2 f(»rt>f ( , 



Note that if a — )■ we get 



t ~ 1 



Vc 



y(t) w v^(t - t ) for t > t 



which shows the branching transportation behaviour. 

Here are the plots of the transportation paths in some cases. 




0.2 0.4 0.6 0.8 1.0 



0.2 0.4 0.6 0.8 1.0 



0.2 0.4 0.6 0.8 1.0 



Figure 4: (a) V{y) = \y\ 2 . (b) V{y) = 8|y|. (c) V{y) = 4|y| - 1 . 
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